A witness for coherent electronic oscillations in ultrafast spectroscopy 
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We report a conceptually straightforward witness that isolates coherent electronic oscillations from 
their vibronic counterparts in nonlinear optical spectra of molecular aggregates: Coherent oscillations 
as a function of waiting time in broadband pump/broadband probe spectra correspond to coherent 
electronic oscillations. Oscillations in individual peaks of 2D electronic spectra do not necessarily 
yield this conclusion. Our witness is simpler to implement than quantum process tomography and 
potentially resolves a long-standing controversy on the character of oscillations in ultrafast spectra 
of photosynthetic light harvesting systems. 



Recently, there has been considerable interest in long- 
lived quantum superpositions of electronic states in pho- 
tosynthetic molecular aggregates and their potential role 
in efficient energy transport in biological conditions |l|,|2j]. 
Evidence for such electronic coherences stems from time 
oscillations in peaks of two-dimensional electronic spec- 
tra (2D-ES) which persist for over 600 fs 0-H]. How- 
ever, coherences between vibronic levels involving a sin- 
gle electronic state exhibit similar signatures in 2D-ES 
[J, lD, |8[ and have been shown to nontrivially affect en- 
ergy transfer [9l4ll|. Although there are additional hints 
that support the interpretation of the oscillations as due 
to electronic states (beating frequencies and comparison 



with all-atom simulations [12|), unambiguous tools to 



experimentally unravel the nature of these oscillations 
are required. A big step has been the observation that, 
under weak coupling to vibrations and negligible coher- 
ence transfer processes, electronic coherences imply oscil- 
lations in off-diagonal peaks of rephasing 2D-ES and in 
diagonal peaks of their non-rephasing counterparts 13], 
whereas general vibronic coherences show up as oscilla- 
tions in any region of either spectra 1A\. However, the 



rephasing 2D-ES of the paradigmatic Fenna-Matthews- 
Olson (FMO) complex exhibits oscillations in both di- 
agonal and off-diagonal peaks, indicating that systems 
of interest may lie in the regime of strong coupling to 
vibrations 15| or exhibit vibronic coherences only [8|. 
Techniques of wavepacket reconstruction [3] or quantum 
process tomography (QPT) [13, LL8( should clearly pro- 
vide an answer at a cost of several experiments. Our 
purpose here is to provide a practical witness for coherent 
electronic oscillations, which is applicable across different 
regimes of weak and strong coupling to vibrations. 

We illustrate the witness by considering the simplest 
molecular exciton model, the coupled dimer 17]. Its 
Hamiltonian is given by H (R) — Tjy+H e i(R), where T/v 
is the nuclear kinetic energy, and H e i (R) is the electronic 
Hamiltonian which depends on the nuclei R, H e i(R) = 
Emn V mn( R )\mn}(mn\ + J{R)(\l0)(0l\ + \01)(10\).\mn) 
denotes the electronic state with m, n excitations in 
the first, second molecules, respectively (m, n £ {0,1}), 
V mn {R) is the corresponding diabatic potential energy 
surface, and J(R) is the coupling between site excita- 



tions. Any pure state |^/) may be expressed in terms 
of vibronic states, that is, product states of the elec- 
tronic (system) and nuclear (bath) degrees of freedom, 
I*) = J2i a i\ e i)\ N i)> for coefficients a it and {|ei),|iVj)} 
electronic and nuclear bases. A reduced electronic de- 
scription of 1^) is obtained by performing a trace over the 
bath, p e i — Tr ?mc (|*)(vl'|). We consider light-matter per- 
turbation in the dipole approximation, H pert (s) = fi ■ 
e(r,s), where fi = Ee=oi,io (v eg \ e }(9\ + M /e |/)(e|)+h.c. 
is the dipole operator, and e(r,s) — ^ PP ,[e p (s — 
ip)e p + c.c.] denotes the pump (P) and probe (P') pulses, 
with e p (s-t p ) = - 7 A= e -^ P (s-h)-(s-t p ) 2 /2a 2 the Gaus _ 



sian time-profile. 



Here, 



A, 



UJ,, 



t v , a, and 



are the 



strength, carrier frequency, center time, width, and po- 
larization of the p-th pulse, respectively. We shall dis- 
cuss PP' spectra Spp>(T) as a function of T = tpi — tp 
(waiting time) [l| , which can be recovered from a 2D-ES 
by integration along both frequency axes (Supplementary 
Material [2l| sec. I, SIT). The main result of this arti- 
cle is: In the Condon approximation and the broadband 
limit (a — > 0), oscillations of Sppi (T) as a function of 
T correspond to coherent electronic oscillations; in this 
limit, Sppi (T) may be expressed solely in terms of re- 
duced electronic states p e i, so oscillations cannot be due 
exclusively to nuclear dynamics. 

The PP' signal may be written as the sum of Sse(T), 
Sesa(T), and Sgsb(T), with separate contributions 
from stimulated emission (SE), excited state absorption 
(ESA), and ground state bleach (GSB) 0. If the initial 
vibrational state is known, each of these terms may be ex- 
pressed as a suitable wavefunction overlap (SI-I 21]). For 
example, let the initial wavefunction (before any pulse) 
be |\T/o(0)) = \g)\v\ 9 ), where \v\ ) is a vibrational eigen- 
state of H vi b, v (R) = T^ + V V {R). Treating the laser 
pulses pertubatively, the first order wavefunction due to 

Pis(ft=l) |*p(a)) = i JZo d3 ' e ~ iHo{s ~ s ' ) {l Jl ■ e p( s ' ~ 
^-p)}|^'o( s ')) j an d the second order wavefunction due to 
both P and P' is |# PP '(s)) = i f^ ds / e - tH « l - s - s ">{^ ■ 
e P '(s')}|* P (s'))- It can be shown that S SE {T) = 



23]) 



(ypp>(s)\g)(g\*pp>(s)) (SI-I [21| and 

Preliminary example. — We will develop some in- 
tuition through an illustration, in which we focus on 
Sse(T). Consider the case where the surfaces of 



the singly-excited diabatic states have the same shape, 
Vio(-R) = Vqi(R) + c, for some constant c (but in general 
V g (R), V f (R) + V e (R) + c for e = 01, 10). It is conve- 
nient to introduce the excitonic basis {\g}, \a), |/3), |/)}, 
which diagonalizes the electronic Hamiltonian at the 
ground state nuclear configuration: H e i(0) = uj g \g}(g\ + 
LJ a \a){a\ + up\0)(fi\ +u) f \f)(f\. Here, \g) = |00) and 
|/) = |11), but in general, \a) and |/3) differ from |01) 
and 1 10) in that they are delocalized due to J(0). Note 
that both | a) and |/3) are coupled in the same way to the 
vibrational bath, and hence they form a decoherence-free 
subspace [25[ . The first order wavefunction "right before" 
the probe pulse may be expanded as |4'p(ii + T)) = 

J2i= a ,0J2i c i,m.(T)\i)\vm)- Since in this case, \i)\vm) 
are eigenstates of the molecular Hamiltonian H (R), 
the excitons are the adiabatic electronic states, there is 
no dissipation in the electronic system, and the values 
\ci,m(T)\ 2 are constants as a function of T, depending 
only on the details of P 28| . The wavefunction "right 
after" the probe at time T is, in the Condon approxi- 
mation, given by, \^ PP , (h +T)) =i J2i= a ,0 £ mn /x mg • 
ep'ip, (ujirn^n)^ 9 ' \vir\ ')ci,m{T)\g)\v£ '), where e p (w) = 
^ e -(w-w p ) a /2 j g £k e p our j er transform of pulse p at 
frequency lo. This expression can be interpreted as a 
wavepacket in the ground state created when the probe 
couples the vibrational levels of the singly-excited states 
to the vibrational levels of the ground state via the elec- 
tric dipole moment, where the amplitudes in the various 
vibrational levels depends on the probe's electric field 
at the given transition energy and the Condon overlap. 
Computing the norm of the resulting wavepacket, 



Sse(T) = Yl (tog ■ e P')(Vjg ■ e P') 

ij=a,/3 

x E (^i^x^W) 

mm'n 
X l P , {uJ im ,gn)e*p> (Ujm',gn)Ci,m (T)c* jm , (T)fl3,) 

which corresponds to sums of interferences between vi- 
brational states of the same and different excitonic states, 
respectively, projecting onto the same vibrational state in 
the ground state. Note that Sse{T) can be written as a 
linear combination of elements of the full vibronic den- 
sity matrix p(T) = \^ P (h +T))(V P (t 1 +T)\. The terms 
(i,m\p(T)\j,m') = c im {T)c* m ,{T) for (i,m) + (j',m') 
correspond to vibronic coherences and oscillate at the 
difference frequency between the |«)|m) and the \j)\m!) 
states. When we consider the broadband (bb) limit of 
Eq. ([IJ , where e(tu) = A for all the w values of interest, 

Sf E (T) = A 2 E (^g-ep^-ep^cUTy^m) 

Crucially, Eq. @ is a linear combination of elements 
of p e i(T) as opposed to the full vibronic space. In fact, 



the terms for i = j correspond to electronic popula- 
tions and, due to the absence of electronic decoherence 
in this example, stay constant with respect to T. The 
term (a\p e i(T)\/3) = J2i c *,»(T)c* f:j (T) corresponds to an 
electronic coherence between \a) and \(3), and shows os- 
cillations at the single frequency ui a p as a function of 
T. Hence, coherent oscillations in S b g E (T) are a wit- 
ness for coherent electronic dynamics. Remarkably, in 
the additional limit where one of the excitons is dark 
(e.g., /j,p g — 0), we have a monomer instead of a dimer, 
and Sly^;(T) is a constant even in the case of large Con- 
don displacements, where there is large vibrational mo- 
tion between pump and probe. This observation for 
the monomer has been previously reported by Yan and 
Mukamel J26|. 

The results above can be interpreted as follows. In the 
Condon approximation, the probe couples only to the 
electronic dipole, so in the broadband limit it acts uni- 
formly across every transition energy, and hence across 
every nuclear configuration within a particular electronic 
state. In general, Sse(T) is a sum of multiple interfer- 
ences among portions of wavepackets at different elec- 
tronic and nuclear configurations. In <Sl^(T), the probe 
opens only two interference pathways (just as in the 
double-slit experiment), via emission from the \a) or the 
|/3) state, insensitive to vibrational dynamics, providing 
a witness for coherent electronic oscillations. 

General case. — The example above readily generalizes 
to include effects of initial thermalized states of the bath, 
ESA and GSB contributions, and non-adiabatic effects. 
In the limit of broadband P (SI-II and III, [H|) and P', 
each of the contributions to S P b P ,(T) are (SI-II, [Jl|). 
Sf E (T) = A 4 E,™> fl4 ' ep/)(/i M • ep)(/x flP • e P )( M . ' 



ebb 



= -A 4 £„„<>« • e P ')(n 



91 

e P')Xijqp( T ), S ESA( T ) ~ ~" l^ijpqWfi 

e p)(v gp ■ e p)(Vjf ■ e P')Xi]q P (T), and 
A4 Ei P (M 9P • e P )ln pg ■ e P )(n gi ■ e P ,){p, ig ■ e P >), where 

(T) = 
Tw{(i|e-^° T (la) (pl®/9 R (0))e^« J | j)} [17LI18L and it 



qg 
S b G b s B (T) = 



Jip\r~gp 

the process matrix x(T) is given by, 
A(Ae- tHoT (\q)(p\® PB(0))e* H ° T \j)} 
is easy to see that S P h P , (T) is invariant under change of 
electronic basis within the singly-exicted states. Here, 
Ab(0) = J2 n Pn\vn ){vn \ is the initial thermal vibra- 
tional ensemble in the ground electronic state. x{T) 
describes the evolution of the electronic system, assum- 
ing that the vibrational system begins in pb(0)- If 
the initial state of the bath can be prepared at p\b(0) 
regardless of the electronic state, as in the impulsive 
limit, an integrated equation of motion can be written 
as Pij(T) = J2i jq pXtjqp( T )Pqp(Q)- As in th e preliminary 
example, S P b P , (T) is a linear combination of entries of 
reduced states p e \ (T) , so oscillations in it are a manifes- 
tation of electronic oscillations, justifying the witness. 

Given an electronic basis, any element Xijqp{T) can 
in principle exhibit oscillations. For a large variety of 
systems, it is however, possible to associate the largest 
amplitude oscillations of x(T) to electronic coherences 



in some basis. In the preliminary example, the lack of 
dissipation implies that Xijqp(T) = S iq Sj P e~' lu ' qpT , so the 
only possible oscillatory contribution to Sppi (T) corre- 
sponds to XaPapiT) = x}i a fi a ( T ) (excitonic coherence). 
In the non-adiabatic case where Vqi(R) ^ Vio{R) + c, 
each electronic state couples differently to the vibrational 
modes. However, in the limit of weak system-bath cou- 
pling, the vibronic states \e)\v ) are still the correct 
eigenstates of H (R) up to zeroth order in the coupling, 
so any oscillations in the signal will still be dominated 
by excitonic coherences. Finally, for intermediate and 
strong system-bath coupling together with a fast bath 
decorrelation timescale, a polaron transformation defines 
an electronic basis {\g), \a), \/3), |/}} that diagonalizes a 
zeroth-order electronic Hamiltonian weakly coupled to a 
renormalized bath ([12( and SI- IV [21(). In this case, the 
highest amplitude oscillations in its S b PP ,{T) would cor- 
respond to electronic coherences Xafs&fsC^) = Xs-s-C^)- 
For more general aggregates, if this were an issue of in- 
terest, a partial QPT could be designed to determine the 
value of specific terms of x{T) [17I. Il8||. 

Numerical examples. — We have performed simulations 
for a monomer, a dimer which exhibits coherent elec- 
tronic oscillations, and an incoherent dimer, where each 
singly-excited site is coupled to a single vibrational mode. 
These three examples illustrate the value of the witness 
(Fig. 1), as all three have oscillatory 2D-ES (Fig. 2), but 
the monomer and incoherent dimer do not have coher- 
ent electronic oscillations. The witness correctly shows 
that only the coherent dimer has a positive witness. The 
simulations include inhomogeneous broadening (ensem- 
bles of 500 molecules with Gaussian site disorder of stan- 
dard deviation 40 cm -1 and, for the dimers, site energy 
correlation 0.8), thermal averaging of initial vibrational 
states according to a Boltzmann distribution at 273 K, 
isotropic averaging, and explicit inclusion of pulses with 
the dynamics. Roughly, there are two energy scales to 
consider, an average coupling J and a reorganization 
energy A, in which case the impulsive limit is set by 
i 3> max(J, A). For these simulations, the pulses are 
within the FWHM=10-20 fs range, and cover the entire 
absorption spectra, respectively (SI-V, 2l|). Fig. 1 shows 



(>Spp' (T))zzzz, the witness averaged at the collinear pulse 
setting zzzz, for about 900 fs (top). We can associate the 
witness oscillations to oscillations of elements in x(T). 
We show a few representative elements of this matrix 
(bottom). Fig. 2 presents snapshots of the rephasing 2D- 
ES, (S(uj T ,T,ujt))zzzz, for a sampling of waiting times T 
between 71.6 and 270.6 fs (left), indicating that vibronic 
coherences manifest as diagonal and cross-peak oscilla- 
tions [29|. Notice that due to strong coupling to vibra- 
tions, the coherent dimer also exhibits oscillations in the 
diagonal peaks, implying the inapplicability of previous 
measures for this case 13|, UA \ ■ As another illustration, 



black is in the right plots. Note that the largest ampli- 
tude oscillations are in the monomer, which cannot have 
coherent electronic oscillations, showing that oscillations 
in peaks in the 2D-ES do not directly translate into co- 
herent electronic dynamics, and hence are not the correct 
witness. 

The witness is positive if, once the dc background is 
subtracted from S PP , (T) , there are oscillations with am- 
plitude proportional to /i 4 , where /i is some estimate of 
an electronic transition dipole moment. If spurious oscil- 
lations due to finite pulse-duration are suspected, a more 
quantitative confirmation is the following: (a) Collect 
traces S PP (T) at several pulse widths er, all roughly in 
the broadband domain, (b) Fourier transform the data: 
S^ P (ui T ) = ijf™dTe iuT Sp b P (T). (c) Locate non-zero 
frequencies of S pp (ujt) corresponding to oscillations be- 
tween discrete states (ignore the dc component). For 
each of these frequencies lot, plot S P b P (uiT) as a function 
of a, and linearly extrapolate to a — > 0. If the obtained 
intercepts are at zero within noise levels, the witness is 
negative . SI-II [21] displays an analytical expression for 
the O(a) correction of S PP (T), providing a theoretical 
basis for this procedure. 

Although the theory has been detailed here for a dimer, 
the witness is applicable to larger aggregates. In the case 
of FMO, due to spectral congestion, it might be fruitful 
to focus on pairs of exciton states at a time, for instance, 
the first and the third exciton states, either via direct PP' 
measurements that cover these transitions exclusively, or 
alternatively, integrating windows of broadband 2D-ES 
corresponding to these two states only, assuming that 
relaxation processes do not occur outside of this spectral 
window. 

We dedicate this letter to the late Bob Silbey who at 
different stages introduced J.Y.Z. to theoretical chem- 
istry and mentored A.A.G in the field. We thank DARPA 
Award No. N66001-10-4060 and EFRC-DOE Award No. 
Ide-sc0001088l 
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Figure 2: (Left) Norm of broadband rephasing 2D-ES for (a) 
monomer, (b) coherent dimer, and (c) incoherent dimer. Di- 
agonal and cross peaks oscillate as a function of time in all 
cases, indicating general vibronic coherences but not necessar- 
ily electronic coherence. Color scale (arbitrary units) is fixed 
throughout. Black dotted circles at left indicate cross peaks 
whose real part amplitude is shown followed in a finer time 
grid in the right panels (varying axes for each sample, arbi- 
trary units). These traces are the wrong witness for coherent 
electronic oscillations. 
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PP' SIGNAL IN TERMS OF WAVEPACKET OVERLAPS 

Consider the situation described in the article, where the total Hamiltonian is given by H — Hq(R) + H pert (s). 
In this section, we will assume Ho(R) to be the same molecular piece as the one described in the article, and 
Hp er t(s) — —fJ, ■ e(r, s) to be the standard light-matter interaction in the dipole approximation, although we consider 
a slightly more general setup, where the electric field is described by three non-collinear beams, e(r, s) = Xln=i[ e p( s — 
t p )e lkp ' r+l ^ p e p + c.c] with different wavevectors k p and phases <fi p . The expressions for Spp>(T) will appear as we 
take the limit of the PE signal to the PP' limit. 

The pulses generate a time-dependent polarization P(r,s) ~ Ti(fip(r)) — J^ fc P(fc; s)e ,fcs ' r on each molecule at 
position r [19(. The allowed wavevectors are the phase-matching directions k — qk\ + rfc2 + sk^ for integers q, r, s, 
and encode different sequences of interactions of the pulses with the molecule. We are interested in the signal S at 
the photon-echo (PE) phase- matched direction kpE = — k\ + &2 + &3, which can be detected by mixing the material 
ensemble emission with a local oscillator (LO) pulse 64(5) travelling along k 4 = kpE, S(r,T,t) = —2^sf_ dse^s — 
£4)64 • P(kpE', t, T, s), where r = £2 — ti (coherence time), T = £3 — £2 (waiting time), and t = £4 — £3 (echo time) [l|. 
Upon repeated collection of S(t, T, £) for many values of time intervals, a 2D-ES can be constructed as a function of T, 
by Fourier transforming the signal with respect to r and t, S(u T ,T, u t ) = J °° dre~ lUTT / °° dte lUtt S{r, T,t) 0, Q. In 
general, oscillations in 5*(t, T, t) and S{io T ,T,u> t ) can be associated to coherent superpositions of vibronic eigenstates 
of H (R), but not necessarily of electronic states [4J. In the article, we paid special attention to the pump-probe (PP) 
limit Spp'(T), which is equivalent to a differential transmission signal, where the first two pulses act as the pump P, 
(ei = €2 = ep, <fii = (f>2 = 4>p), the last two as the probe P' (£3 = £4 = ep> and 03 = 04 = <fip>), t — t — 0, and P and 
P' are well separated (i.e., T » a), Spp>(T) can be recovered from the 2D-ES as an inverse Fourier transform at zero 
frequencies, Spp>(T) = 5ps(0,T, 0) = tAt /^ dw T ^^ (lwti?(w T , T, Wt). This limit justifies the form of H pert {s) 
given in the article, which only consists of two pulses. 

The starting point is the expression for Sppi (T) , 

/oo 
ds'el{s' - t 4 )e 4 ■ P(k PE ; 0, T, 0). (SI) 

-00 

We shall derive a wavepacket overlap formula for Spp>(T) assuming that P and P' are well separated, T ^> er, 
analogously to the doorway- window approach [l|. First, we conveniently define the following wavefunctions: 

|tto(*)> = e-^-^lMs')), (S2) 

/•OO 

|*p(a)>=t/ ds / e- iH ^-^{ f ,-ep(ep(s / -tp)+c.c.)}\^ Q ^)), (S3) 



\*pp>(s)) =1 ds'e-^ s -^{ii-ep,(ep,(s'-tp,)+ C . C .)}\^p,{ S ')), (S4) 

J — OO 

POO p8 

\Vpp(s)) = {i) 2 / ds' / ds' / e- tH ^ s - s '){fi-ep(e P (s , -tp) + c.c.)} (S5) 



— OO J — CO 



xe -iH (s'-s") {fJi . ep(£p(s » _ tp) + c.C.)}|*0(«')), 

/oo ps' 

ds' ds"e- lH ^ s ' s '\^-ep,{€ pl {s' -t P ,)+ c.c.)} (S6) 

-OO J — CO 

xe- fl » (s '- s 'V ■ e P ,{ep,{s" - t P >) + c.c.)}|*o(*")>» 

/OO /*s' 

ds 1 ds"e- lH ^-^{^-e P ,(ep,{s' -t P ,)+ c.c.)} (S7) 

-OO J —OO 

xe -iH {s'-s") {yL . ep ,( ep ,( s " _ fp/ ) + c.C.)}\* P p(s")), 



which are valid for s ^> tpi (after the envelopes of the pulses have considerably decayed), and which correspond to 
the processes indicated by their subscripts, i.e., \^fppp>p'{s)) corresponds to the fourth order wavefunction (0(A)) 
resulting from two actions of P and two of P' (see Fig. SI). 

Eqs. (|S2]) - (fS7)) allow for a calculation of P(k PE \ 0, T, 0) and hence of S PP > (T) via Eq. (jST) . Note that, as opposed 
to a general PE signal, Spp> (T) does not depend on the phases of the pulses because </>i — <p2 = <Pp an d <fe = (f>4 = 4>p> ■ 
The phase-matching condition kpE = — fei + &2 + k% together with the rotating-wave approximation indicate that 
for each wavevector +kj(—kj), the pulse j acts with the term £j(e*), exciting (de-exciting) the ket or de-exciting 
(exciting) the bra. Collecting all the terms result in S PP >(T) = S S e(T) + S E sa(T) + S GS b(T): 

/CO 
dt'^ Ppi {t')\{e pl {t' - t P ,)ti ■ e P ,}|*p(t')> 
OO 

/OO ft' 

dt' / ds'(Vp(s')\{-iep,(s' - t P >)ix ■ ep/}e iHo( *'-*' ) {(*)cJ«(t' - tp>)pt ■ e pi }\^ P {t')) 
-oo J —co 

= 2{$pp,( S )\g)(g\$ P p,( S )), (S8) 

/CO 
dt'^> P {t')\{€ P ,{t' - t P ,)n ■ ep/}|*p P /(0> 
-OO 

/OO ft 

dt' / ds' {^ P {t')\{{-i)e PI {t' - t P .)fj. ■ e P }e- lHo(t '- s,) {ie P (s' - t P )ft ■ e P }\* P (s')) 
-OO J — OO 

= -2<*PP/(s)|/)</|*pp/(*)) ) (S9) 

/OO 
dt'{(y PP {t')\{e P ,{t' - ti»)i* ■ e P >}\*p>tf)) + (*o(*')|{cp/(t' - tp>)p ■ ep>}\^ppp>{t'))} 
OO 

/CO 
dt'{(*pp(t')\{(i)ep,(t' - t P i)n ■ e P ,}\*p,(t')) + (* (t')\{(i)e P ,(t' - t P >)n ■ e P ,}\* PP p,(t'))} 
-co 

= 2%t{{*pp{s)\g){g\^p lp ,{8)) + (V (s)\g){g\*PP P >P>(s))}. (S10) 

where again, s 3> tp', and otherwise, the signals are independent of s. This can be easily understood in physical 
terms: After the action of the pulses, the wavefunctions still evolve according to Hq(R), but the overlaps do not 
change in time. Eqs. (|S8|) — (IS10[) are in the spirit of wavepacket approaches to PP' spectroscopy ja-llll]. 

In order to gain additional insight, we interpret the formulas in terms of differential transmission by enumerating 
all the possible absorption and emission processes which are quadratic in P and P'. P promotes a wavepacket from \g) 
to |\I / p(s)), a superposition of wavepackets in \a) and \j3). P' acts on this state, creating |^pp/(s)), a superposition of 
wavepackets in \g) and |/). Naturally, the photons emitted in SE correspond to the squared amplitude of (g\^ppi(t)), 
whereas the ones absorbed in ESA are associated with the squared amplitude of (f\^>pp>(s)), hence providing an 
intuition for the expressions for Sse{T) and Sesa(T). Finally, Sgsb(T) can be thought as accounting for the 
"leftover" SE processes, namely, overlaps between wavepackets created by pulses at different times. After P and P' 
, the total ground state wavepacket is (g|*(i)) = (g\(\^ a (t)) + |*pp'(*)) + I*ppW) + l*P'P'W) + I*ppp'P'W) + 
higher order contributions). Collecting wavepacket overlaps which are quadratic in both pulses yields Sse(T) + 
Sgsb (T) . "Leftover" ESA processes do not contribute to Spp> (T) because they do not fullfill the PE phase-matching 
condition (they appear in double-quantum coherence spectroscopy, for instance). 

Thermal effects follow from averaging the signals corresponding to initial states \^o(tp)) sampled according to a 
Boltzmann distribution. 

GENERAL EXPRESSIONS FOR S PP ,{T) IN VIBRONIC BASIS 

In order to manipulate the wavepacket overlap expressions, it is convenient to define two vibronic bases: 
• The vibronic eigenbasis of Ho(R),{\g, Vn ), |C), l/> v n )}, which satisfy 

H (R)\g, vi 9) ) = u> gn \0, (Sll) 

H (R)\0 = w c |C), (S12) 

H {R)\f,vi f) ) = u /n |/,i#>). (S13) 
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Figure SI: Feynman diagrams for the wavefunctions denned in Eqs. (|S2[) - (|S7 



where |£) corresponds to the singly excited state manifold. 

• The tensor product basis {\m, i/„ )}, where {|?n)} denotes electronic states in an arbitrary electronic basis (for 
instance, the excitonic one), and {\vn )} refers to vibrational eigenstates of the ground vibrational Hamiltonian, 

H V ib,oo(R) — Tjy + Voo(R). 

Note that we can always write states in the vibronic eigenbasis in terms of the second one: \g, Vn ) stays the same, 

10 =E m „<m,^ ) |C)K4 ff) >, and \fA?) = EJ/,*# ) X/>*# ) l/,»# ) > = Y m ($ ) \v i #)\f,v$)- 

Using both bases, the process matrix x(T) affords a compact representation, 



Xi JqP (T) = X lJ9P (r)-Tr„ uc {(*|e-^ oT (|g)(p|®PB(0))e^ T |j)} 



(S14) 



££'nn' 



Our goal is to write Spp>(T) for arbitrary bandwidth in a similar style, so that in the broadband limit, we can 
identify it as a linear combinations of elements of x(T), hence proving Eqs. (|S22[) — (IS24|) in the article. We start by 
rewriting Eqs. (lS2J) - (|S7f in the vibronic bases: 
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(S15) 



(S16) 



|3>I*# ; >> (S17) 



X ( 1 - Erf ( i<T ^ UJL + "V'.c) + i^L - ^C'.gn') 
x (l-Erf( t ' T(K+ ^^ +( ^'^" ) ^ 1.9)1^). (S18) 

where the Erf functions appear due to pulse overlap. Eqs. ([58]) — (|S10|) together with Eqs. (|S15[) — (|S18[) yield: 

Sse(T) = ^2(fJ> gi -ep,)(n qg ■ e P )(n gp ■ e P )(n jg ■ e P >) ^ tp>(ugn>£)ep(wc,gn)ep(ugn,<;')ep'(v<;',gn>) 

ijqp CC'nn' 

xpne-^-^'^iiMjHoickA^A'Hc'KcU^), (si9) 

Sesa(T) = -^2(n fi -ep>)(fJ, qg -ep)(iJ, gp -ep)(iJ, j fep>) 

ijqp 

X H (^"Vm^^mVi' >eP'(w/m,c)ep(wC,9n) ? p( w 9n,C')eP'(wC',/m) 

££'nn / n // m 

x^e-^-^^^^icxck.^Xp.^lC'XC'li,^), (S20) 

^Gss(r) = 23t^2(ti gi -e P ,)((j. qg -ep'Xflgj, ■ e P )(n jg ■ e P ) l-j ^ ep'(uj gn 'x)ep'(uj Cgn )ep(uj gnX ')ep(uj Ctgn/ ) 

ijgp CC nn ' 
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+ Pll ( i ,^ ) |C)(Cl9,^ ) >((p,^ ) IC')(C'|j,^ ) )e^"'- T )* 

f ia{{-UJ gn . x - u L ) + (u} Cgn - u L )) \ \ / Erf Acr((-ay, c - u L ) + (^ C ',g«' ~ ^l)) 
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The expressions above can be intuitively understood in terms of the double-sided Feynman diagrams in Fig. S2. The 
expression for GSB consists of a sum of terms corresponding to two types of Feynman pathways which are different 
in general. 

In the broadband limit where ep(cj) = epi(u>) = A, many sums above collapse through resolutions of the identity, 
and we straightforwardly recover the expressions in the article, 

Sf E (T) = A 4 J2(n gi ■ e P ,)(n qg ■ e P ) (S22) 

ijpq 

X (Vgp ' e P)(Mj9 ' e P')Xijqp(T), 

S b E b s A (T) = -X^fafi-ep'X^-ep) (S23) 

ijpq 

x (m 3P ' e p)(Mj/ ' e P >)xijqp(T), 
S b G b s B (T) = A 4 £> OT • e P )(» pg ■ e P ) (S24) 

ip 

x(l*gi-ep>)(Hi g -ep>). 

EIn this limit, as highlighted by the T- independent form of Eq. (|S24[) . the two types of GSB pathways yield the same 
stationary background to the signal (caused by copies of the initial stationary wavepackets in the ground electronic 
surface) . 

In the practical case where the pulses are broad, but not infinitely sharp in time, we can expand, Sppi (T) = 
S P b P ,(T) + S { PP ,(T) where SP P ,(T) = s£l B (T) corresponds to corrections of O(a), which originate from the Erf 
functions in the GSB signal: 
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Figure S2: Double-sided Feynman diagrams for general PP' signal. 
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SE and ESA processes only contribute to corrections of 0(a 2 ) via the Gaussian spectral profile of the pulses. 



REQUIREMENT OF BROADBAND PUMP P 

Although the conclusions of the premilinary example in the article hold even in the case of narrowband P', we also 
require broad bandwidth for P for two reasons: 



1. N on- stationary GSB contributions. Eqs. (|S5[) and (|S7J1 show that in the limit of broadband P, this pulse promotes 
a wavepacket to the excited states and immediately back down to \g), yielding a wavefunction |\&pp(£)) that 
is proportional to the original Y$o{tp)) before any pulse (also see Eq. IS24|) . In this limit, as emphasized in 
the previous section, Sgsb(T) is a constant background as a function of T, giving the opportunity to identify 
Spp'(T) as a probe for singly-excited state dynamics. Under a narrowband P, this no longer holds, as shown 
by Eq. (|S21[) . which depends on T in general. In this case, |\l/pp(i)) will be a non-stationary wavepacket in the 
ground electronic surface, which will manifest as time-evolving overlaps both in ( x I'pp(s)|(7)(g| , I'p<p<(.s)) and in 
(*o(8)\g){g\*ppp> P >(8)) (seeEq. (SED). 



Consistency with QPT. As mentioned in the article, the initial states prepared under a broadband P are of the 
form p(0) = \q)(p\ ®ps(0), where pb(0) = J2 n Pn\ v ™ )( v n | is the initial thermal ensemble of vibrations in the 
ground electronic surface for all |<7)(p|. Under a narrowband P, it is not possible to prepare initial tensor product 
states between the system and a fixed bath pb, so Spp>(T) can no longer be written in terms of elements of a 
single x(T), and the equation p(T) = x(T)p(0) loses its meaning. 



POLARON TRANSFORMATION 

Here, we summarize the essential features of the polaron transformation used in the arguments of the article. We 



closely follow the works of Silbey, Harris, and coworkers [12H14I ]. Consider the approximation where the diabatic 



2 r>2 

potential energy surfaces are given by harmonic wells along each nuclear coordinate, V\q{R) = Eiq + J2 n m "^" " + 
y / '2m n ujf l gio, n R n , and Vqi(R) has the same form except for the substitution 10 — > 01 in the subscripts, whereas V g (R) 
and Vf(R) have arbitrary shapes. Here, w„ denotes the n-th mode frequency, whereas the displacements <?io,n denote 
linear couplings of the electronic system to the nuclear bath. Define the harmonic oscillator creation and anhilation 

operators in the usual way bk — \f^f-x n =F / 2 m = af~' an< ^ a ^ so ^e generator G = Y^ n (Pn ~ &n)(<7io,n|10}(10| + 
5oi,n|01)(01|) such that U — e G corresponds to a unitary transformation of the full-polaron transformation [l4j . 
It follows that H(R) = e G H (R)e^ G = H (R) + Hi(R), where H (R) = T N + H e i(R) is our new zeroth-order 
Hamiltonian, and Hi(R) is the perturbation term, whenever it is small compared to Hq(R). To make a connection 
with the previous notation, we explicitly write H e i = J2mn V mn \mn)(mn\ + J(|10)(01| + |01)(10|), where, 

Vw - ti w - ^uj n g Wn + 2_^ ^ ' 

n n 

V — TP S^ , , «2 , V* m nU n R n 

Vai - ti 01 - 2_^u n g 01n + 2_^ ^ , 

n n 

J = J{w), 

w = exp I ^(flfio.n ~ 9oi, n )(bi -b n )\ , 

(w) = Tr(wp B {0)) 

= ex P I - Z^ coth 2 J ' 

The expressions /3, p B (0) = J]„ J2 r — z n \l)' T ""^ > Zn ^) = iimE^Ti) ' ' w ™' r ^ label the inverse temper- 
ature, the initial thermal ensemble of vibrations, the partition function of the n-th oscillator, and the r-th eigenstate 
of the n-th harmonic oscillator, respectively. J can be interpreted as a renormalized site-site coupling due to phonon- 
dressing. Furthermore, Hi(R) = J(w — (iu))|10)(01| + h.c. Jang advises to consider the smallness of the quantity 
Jyl — w 2 as the figure of merit for the validity of perturbation theory, and hence for the usefulness of the polaron 
transformation. Cao and coworkers note that the accuracy of the polaron transformation is guaranteed only in the 
scenario of fast bath decorrelation compared to the other relevant timescales [15| . 

If the dynamics of all the degrees of freedom are governed by H${R) alone, the electronic system is effectively 
uncoupled from the nuclear bath. The diagonalization of Hq(0) yields polaronic states {\g), \a), \/3), |/)} which satisfy 
Xijqp(T) — 5i q Sj P e~' lu ' qpT . As can be easily checked, Eqs. (|S22| — [S24J1 are invariant under change of basis. Hence, if 
H\(R) can be guaranteed to be a small perturbation for Hq(R), to zeroth-order in Hi(R), the coherent electronic 
oscillations in S P b P , (T) correspond to electronic coherences in the polaronic basis. 

The steps above have been outlined for the full-polaron transformation, but the conclusion can be easily seen to 
hold whenever the total Hamiltonian can be repartitioned into a large contribution and a small system-bath coupling. 



Examples include the variational polaron transformation [12LI13J. which interpolates between weak and strong coupling 



between the original system and bath, as well transformations which include anharmonicities in the diabatic potential 
energy surfaces (quadratic coupling between the original system and bath, |16|). 
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Figure S3: Inhomogeneously broadened absorption spectra (solid red) with pulse spectral profiles |e p (u>)| on top (dotted black). 

DETAILS OF THE NUMERICAL SIMULATIONS IN THE ARTICLE 

We have performed computational simulations 20] for absorption spectra, PP' signal, and rephasing 2D-ES for a 
monomer, a dimer which exhibits electronic coherent oscillations, and an incoherent dimer. For their Hamiltonians, 

1 1 T 1 r 1 1 T 7 t \ 7-1 U -rix — ^-mn x ) OJ (y — A^n y ) . 

we choose harmonic diabatic surfaces parametrized by V mn (x,y) = h mn -\ ! ^ : 1 — ^ — 2 — > where 

x and y are scaled nuclear coordinates, E mn are site energies, u mrl x i y \ are oscillator frequencies and A mn x r y \ are 
electron- nuclear couplings [6|-|8j, |l7J, [l8| . The parameters for the calculations are listed in Table 1 . We assumed that 
Vf(x, y) = Vxo(x, y) + Vqi (x, y), and that the carrier frequency of all the pulses is lo p — u)l. The dimers are such that 
the dipoles are oriented 90 degrees from each other, and the ratio between their norms is 1:3. 

Fig. SI shows inhomogeneously and rotationally averaged absorption spectra (solid red) for the examples used in 
the numerical simulations of the article. The pulse profiles (dotted black lines) which roughly cover the strongest 
vibronic transitions with similar amplitude, giving a qualitative idea of what broadband means in practice. In these 
examples, FWHM = 2y / 2ln2a = 10, 18.7, and 12.5 fs, for the monomer, coherent dimer, and incoherent dimer cases, 
respectively. 



TABLE 1. Parameters of Simula 
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£00 /cm- 1 











(£10 -uj l ) /cm- 1 


-125 


-300 
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(E 01 -WiJ/cm" 1 
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J/cmT 1 
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100 


10 


Uqq.x = u>oo,y/cm 


100 
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ww,x = woi.a/cm -1 
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wio : y = woi.y/cm^ 1 
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Aoo.x/cm 1 / 2 
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A^^/cm 1 / 2 
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10 
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